DESY 05-002 
SFB/CPP-05-02 
January 2005 



hep-ph/0501178 



Mellin Moments of the Next-to-next-to Leading Order 
Coefficient Functions for the Drell-Yan Process 
and Hadronic Higgs-Boson Production 



Johannes Bliimlein a and Vajravelu Ravindran ' 6 



a Deutsches Elektronen- Synchrotron, DESY, 
Platanenallee 6, D- 15738 Zeuthen, Germany 

b Harish- Chandra Research Institute, Chhatnag Road, 
Jhunsi, Allahabad, India. 



Abstract 

We calculate the Mellin moments of the next-to-next-to leading order coefficient functions 
for the Drell-Yan and Higgs production cross sections. The results can be expressed in 
terms of multiple finite harmonic sums of maximal weight w = 4. Using algebraic and 
structural relations between harmonic sums one finds that besides the single harmonic 
sums only five basic sums and their derivatives w.r.t. the summation index contribute. 
This representation reduces the large complexity being present in x-space calculations and 
is well suited for fast numerical implementations. 



1 Introduction 



The Principle of Simplicity 1 is one of the guiding principles in physics [1]. Whenever possible 
one seeks for as simple as possible expressions, not only to obtain the result in a more compact 
form, but also to reveal the basic structures behind. This applies also to complex computations 
in particle physics. Without achieving suitable simplifications it is often impossible to undertake 
even more involved calculations, such as one order higher in the coupling constant in perturbation 
theory. 

In the present paper we seek for a simplification of the 2-loop coefficient functions for the 
Drell-Yan (DY) process both for unpolarized and polarized nucleons, and associated to it, to 
those for hadronic scalar and pseudo-scalar Higgs boson production in the heavy-mass limit. 2 
The pioneering 2-loop calculations 3 in this field [5,6] were carried out during the last 15 years 
and paved the way to understand single-scale quantities in Quantum Chromodynamics (QCD), 
beyond the level of the pole terms. They form one of the milestone in the history of QCD. 
The calculations were performed in x-space using the QCD-improved parton model. Here x 
denotes the fraction of the momentum of a radiated particle to that of the source particle. For 
the representation of the Wilson coefficients a set of up to 77 functions is needed, not counting 
those, which lead structurally to the same Mellin transform [7]. 

With the present collider Tevatron [8] and the upcoming large hadron collider, LHC [9], 
at CERN, the need for precise predictions from theory has become more and more important. 
These experiments are aimed at not only to discover new particles but also to measure various 
parameters of the Standard Model at higher precision to confirm the predictions of the theory. 
The Higgs boson of the Standard Model has yet to be discovered at these machines. Needless to 
mention that one also hopes to discover various new particles predicted in scenarios beyond the 
Standard Model, as for example supersymmetric particles or leptoquarks. The measurement of 
the masses and couplings of newly discovered particles likewise the known particles require the 
precise knowledge of the QCD corrections. 

On the theoretical side the relevant observables have to be calculated at high accuracy. The 
QCD corrections have a considerable impact both in the discovery channels as well as for the 
precision measurements. Theoretical uncertainties emerge from ultraviolet (UV) renormalization 
and mass factorization scale dependence. These uncertainties can only be diminished extending 
the calculations to higher orders. 

In the differential and total production cross sections for the Drell-Yan process as well as 
for Higgs production at hadron colliders the UV-renormalization scale dependence is due to the 
strong coupling constant. On the other hand, processes like Higgs production are associated 
to heavy quark final states such as tt, bb, or inclusive heavy quark productions. They receive a 
UV-renormalization scale dependence also from heavy quark mass effects. 

A second source of uncertainty is due to mass factorization. Both the partonic cross sections 
and the non-perturbative parton densities depend on the factorization scale. Since the partonic 
cross sections are computable order by order in perturbation theory, their dependence on the 
factorization scale is known completely. Similarly, one can fully determine the factorization 
scale dependence of the parton densities through solving the associated renormalization group 
equations (RGE) which govern their evolution. To diminish both effects significantly usually 

1 Pluralitas non est ponenda sine neccesitate. 

2 For the calculation of the Higgs and pseudo-scalar Higgs production cross section including mass effects to 
NLO see [2] and the recent reviews [3]. 

3 For NLO results in the case of Higgs-production see [4]. 
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next-to-next-to leading order (NNLO) QCD-corrections have to be performed. 

In recent years, there have been several significant developments in achieving very precise 
theoretical results to match with the accuracy of the experiments. With the advent of the 
NNLO non-singlet and singlet anomalous dimensions [10] NNLO predictions for the Drell-Yan [5] 
and Higgs-pro duct ion cross sections [6] are possible. To get some feeling for the theoretical 
uncertainty coming from the scale-dependence we consider the following two processes: the 
total Higgs production cross section and the cross section for Higgs production associated with 
a top anti-top quark pair. The total Higgs production cross section is one of the most important 
processes at LHC from the discovery point of view. This production cross section receives 
dominant contributions from gluon-gluon fusion through the top quark loop. Since the leading 
order correction is of 0(a 2 s ) and also involves two gluon densities in the initial state, the NLO 
correction is inevitable to reduce the scale uncertainty. Unfortunately the NLO correction is 
about 80 — 100% larger than the LO correction, casting doubt on the reliability of perturbation 
theory. Interestingly, the NNLO corrections to the process not only give smaller contributions but 
also reduce the scale uncertainty significantly, which improves the stability of the perturbative 
result. Similarly, Higgs production associated with a top anti-top pair suffers from huge scale 
uncertainties of the order of 100 — 200% at the LO level. The NLO correction stabilizes the 
result significantly. The Drell-Yan cross section, which is known upto NNLO [5], not only tests 
the reliability of perturbative QCD but also reduces the uncertainties coming from theory in 
order to make background studies more reliable for new particle searches and physics beyond 
the Standard Model. 

The NLO and NNLO results are technically complicated and they result in large expressions 
involving a large number of functions. For example the perturbativly computable coefficient 
functions at NNLO for the Drell-Yan or Higgs total cross section depend on nearly 80 individual 
functions. These belong to the class of Nielsen integrals [11] with a variety of partly complicated 
arguments and products with logarithms or the weights 1/(1 ±x). For the physical cross sections 
these are convoluted with the respective evolved NNLO parton densities. Since the different sets 
of parton density are usually available in form of numerical codes, the Mellin convolution is done 
performing the respective integrals numerically. Because of the complexity of the expressions 
involved, the numerical computations are usually slow. 

We will discuss here an alternative method which is faster to compute the cross sections. 
It involves the use of the Mellin transformation technique both for the Wilson coefficients and 
the parton densities. The evolution equations of the parton densities can be solved analytically 
in Mellin-space. If the respective Wilson coefficients are also known in Mellin-space one may 
calculate the cross sections by a single numerical integral around the singularities of the product 
of both in the complex plane. This requires the Mellin-representation of the Wilson coefficients 
in the variable N and their analytic continuation to complex values of N [12]. 4 For integer values 
of N the Mellin transformation of the Wilson coefficients leads to finite harmonic sums [7, 14]. 
We use various algebraic relations between these sums [7, 15-17] to simplify the expressions. 
Furthermore, structural relations between the Mellin transforms [18] are used to reduce the 
number of basic functions further. The use of finite harmonic leads to a synchronization of the 
expressions. Several complicated harmonic sums present for individual contributions disappear 
in the final result. A brief summary of the results presented below was given in [19]. Similar 
investigations as performed in the present paper are carried out for the Wilson coefficients for 
deeply inelastic scattering off unpolarized and polarized targets [20] and the unpolarized and 

4 The analytic continuations of the basic Mellin transforms which contribute to the 3-loop anomalous dimen- 
sions have been given in [13] recently. 
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polarized time-like fragmentation functions [21] to NNLO. 

The paper is organized as follows. In section 2 we summarize aspects of the calculation of 2- 
loop Wilson coefficients in x-space. The Mellin moments are discussed in section 3. In section 4 
a brief survey is given on multiple harmonic sums and their relation to Mellin transforms. The 
algebraic relations between the finite harmonic sums used in the present paper is discussed in 
section 5. In section 6 the basic functions representing the 2-loop Wilson coefficients for the pro- 
cesses considered in the present paper are identified. In the appendices A. 1,2 the explicit results 
are summarized for the Mellin transforms of the NNLO Wilson coefficients for the unpolarized 
and polarized Drell-Yan process and appendix A. 3,4 contains the corresponding expressions for 
the hadronic scalar and pseudo-scalar Higgs boson production processes. 

2 Coefficient Functions 

Due to mass factorization, hadronic cross sections such as those for the Drell-Yan process [5] 
and hadronic Higgs boson production [6] can be expressed in terms of Mellin convolutions of 
the perturbativly computable coefficient functions A ab (x, Q 2 , /i 2 ) and non-perturbative parton 
distributions, f a (x,fi 2 ), of incoming hadrons, 



Here the sum over a,b = q,q, g is implied. Q 2 is the mass squared of the Drell-Yan pair or the 
Higgs boson. The parameter x is a scaling variable defined by 



where Pi, P2 are the momenta of incoming hadrons. We have set both the UV- and factorization 
scales to be equal to /z 2 . Though the parton densities are not calculable in the perturbative QCD, 
their evolution with respect to the factorization scale /i 2 is computable using the renormalization 
group equations (RGE). Hence higher order corrections to hadronic reactions enter through two 
sources viz, the coefficient functions and the RG equations of the parton distribution functions. 

The coefficient functions are computed from the partonic cross sections in powers of the 
strong coupling constant a s (fj, 2 ) in perturbative QCD. They are expressible in terms of the 
scaling variable z = Q 2 / s, where s is the center of mass of incoming partonic system. Q 2 is a 
large invariant mass characteristic for the process. The lowest order contributions are of 0(a°) 
and the Wilson coefficients are obtained from the Born-diagrams. In the case of the Drell- 
Yan process the lowest order contributions are due to quark anti-quark annihilation, likewise 
gluon-gluon fusion in the case of Higgs boson production through the top-quark loop. Hence 
the results are just proportional to 5(1 — z). To the next-to-leading order (NLO) one encounters 
virtual corrections as well as real gluon emissions to the Born-processes. In addition, gluon- 
initiated processes also contribute in case of the Drell-Yan process as quark-initiated processes 
in case of Higgs boson production. The virtual processes are in general UV-divergent due 
to the loop-corrections. Since one is dealing with light partons inside the incoming hadrons, 
the real emission processes suffer from collinear divergences. In addition, the standard soft 
divergences appear in both real emission as well as virtual contributions to Born processes. 
All these divergences are regularized using n-dimensional regularization. The UV-divergences 
are removed by standard UV-renormalization. This introduces the renormalization scale fi 2 R 




5(x — zxix 2 ) ■ 



(1) 




S=(P 1 + P 2 ) 2 , 



(2) 
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both through a s (fi 2 R ) and \n(Q 2 / /i^-terms. The soft divergences cancel when virtual and real 
emission contributions are added. The remaining collinear divergences are removed by the mass 
factorization procedure. This introduces a new scale called factorization scale jj? F . Finally one 
ends up with the regular functions as ln(x), ln(l — x), polynomials in x and distribution functions 
such as (ln m (l - x)/{\ - x)) + with m — 0, 1; 5(1 - x). 

Beyond NLO, the computation becomes more complicated. First of all there are many more 
processes that contribute to this order. For example one finds, two-loop virtual contributions 
to the Born terms, one loop virtual corrections to the NLO contributions, double real emissions 
to the Born contributions, and single real emissions to the NLO processes. In addition several 
new type of processes start contributing from NNLO level onwards. Due to this, the number 
of integrations involved increases enormously and new functions emerge. The two-loop virtual 
corrections to the Born process generate only constants such as £ 2 , C3? C±- But one-loop correc- 
tions to NLO processes and real emissions to Born-, NLO- and the new NNLO-processes lead 
to hypergeometric functions F 2)1 which upon integration generate further higher functions. The 
hypergeometric functions usually result from the one-loop box to NLO processes and also from 
angular integration of the real emission processes at the NNLO level. Higher functions are also 
generated when NLO-splitting functions are convoluted with the Born- and NLO-cross sections. 
Such convolutions are required in order to perform mass factorization. The NNLO results are 
finally expressible in terms of polynomials in x, logarithms and Nielsen integrals S njP (x) [11], 
both with various argument-functions, defined by 

The weight w of these functions is defined by w = p + n, where any power of a logarithm counts 
for w = 1 as the case for the denominators 1/x, 1/(1 +x) and 1/(1— x). In a product of functions 
the weights of the factors add. The standard polylogarithms are related to the Nielsen integrals 
by: 

t • / \ ^Li n+1 (a;) 

U n (x) = = S n . 1A (x) 

rid* 



PL ft y 

/ _ ln n - 2 (^)ln(l -zx) for n>2. (4) 
Jo z 



(n-2) 

One relates the logarithms to the polylogarithms by simple differentiation: 

dLi 2 (±x) 



d\n(x) 



Lii(zbr) = - ln(l =p x) , 



Lio(x) = -4-. (5) 
Similarly, one reduces the weight of a Nielsen integral by differentiating w.r.t. \a(x) : 

In addition, distributions as 6(1 — x) and (ln m (l — x)/(l — x)) + with m — 0, 1, 2, 3 contribute. 
These functions and their Mellin transforms can be found in Ref. [7]. 
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The hadronic scattering cross sections are obtained using Eq. (1). One performs the inte- 
gration over Xi and x 2 after folding the perturbativly computed coefficient functions with the 
appropriate parton distributions. This may involve further evaluation of various Nielsen inte- 
grals and increases the complexity of the numerical evaluation of the hadronic cross sections. In 
the next section we will study the structure of these corrections in Mellin space using algebraic 
identities which relate the resulting finite harmonic sums. We will present an alternative treat- 
ment of the evaluation of the total cross sections upto NNLO by working in Mellin space. Such 
techniques have been used in the past to compute deep-inelastic scattering cross sections. They 
are also found to be most suitable for various resummation programs [22]. 5 

3 Mellin Moments 

The Mellin-transform [24] of a given function F(x) is defined by 

M[F](N) = f 1 dxx N - l F(x) . (7) 
Jo 

The Mellin convolution of two functions F 1 (x),F 2 (x) is given by 

[Fi <g> F 2 ] (x) = f 1 cba f 1 dx 2 F 1 (x l )F 2 (x 2 )5{x - x ± x 2 ) . (8) 
Jo Jo 

M [[Fi <g) F 2 ] (x)](N) reduces to the product of Mellin moments of F 1 (x) and F 2 (x), i.e. 

M[Fi <g> F 2 }(N) = M[Fi](JV) • M[F 2 ](N) . (9) 
One may generalize this property for the convolution of m functions : 

[Fi® F 2 ® ... ® F rn ](x) = [ dxi [ dx 2 ... [ dx m [ dzF 1 (xi)F 2 (x 2 )...F m (x m ) 

Jo Jo Jo Jo 

x 5(x - xix 2 ...x m ) , (10) 

m 

M[Fi <g) F 2 ® ...F m ](N) = J] M[F k ](N) . (11) 

k=i 

Indeed, already the multiple convolution of rather simple functions in x-space may lead to com- 
plicated expressions, cf. Ref. [25]. Contrary to that, the representation (11) is straightforward. 

Due to mass-factorization the QCD-improved collinear parton model relates the hadronic 
cross sections via Mellin convolutions to the partonic cross sections and the parton distribution 
functions. Hence the cross section in the Mellin-iV space becomes 

M[a](N,Q 2 ) = M[f a ](N,/i 2 )M[f b ](N,v 2 )M[A^N,^ . (12) 

The Mellin moments of these functions can be analytically continued [12] to complex values of 
iV so that one can use various analyticity properties of these functions in complex iV-space to 
evaluate them efficiently. To retrieve back the full x— dependent result for (12), we have to take 
the inverse Mellin transform as a numerical contour integral around all singularities in N. 

5 Heavy flavor contributions in Mellin space were treated in [23]. 
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For our analysis, the starting point is Eq. (8) with given parton densities f a (x, fi 2 ) and known 
coefficient functions A ab (x,Q 2 ) computed upto NNLO in perturbative QCD. We then compute 
the Mellin moments of these functions in iV-space and analytically continue them to complex 
iV-space. At the end, we use Eq. (12) and perform the inverse Mellin transformation to arrive at 
the results in x-space using a suitable contour in the complex iV-space. Since only one integral 
has to be carried out numerically, the evaluation can be performed very fast. 

Before we study the Mellin moment of the coefficient functions, we would like to make a few 
remarks on the parton densities. As is well known, the parton densities are fitted as functions of 
x using the available deeply inelastic scattering data. The parton distributions are determined 
at some scale Ql along with the QCD-scale, Aq CD . The renormalization group equations for 
mass factorization in x-space 

» — -j-2 = 4?r 2^ J x —Pab{z,a s (n))f b/P (z,n) a = q,q,g (13) 



b=q,q,g 



relate the parton densities at Qq to those at the scale fj 2 . The splitting functions P ab (z, a s (/j, 2 )) 
are computable order by order in perturbation theory: 



P ab (z,a s (» 2 )) = £ P%\z) . (14) 

Instead solving the integro-differential equations (13) one may Mellin-transform these equations 
to 

2 dM[f a/P ](N,fi 2 ) a s (fi 2 ) ^ 2^ A/r[ . i/» r 2^ /-.^ 

» j"2 = — A 1^ M[P ab }{N,a s {fi ))M[f b/ p\(N,fi ) . (15) 

W 471 " b=q,q,g 

The solution of (15) for M[/ a /p](JV, fj 2 ) is straightforward as it is now just a first order differential 
equation. The parton densities in a;— space are obtained from the solutions M[/ a /p](JV, fi 2 ) by 
an inverse Mellin transformation. Since we are dealing with cross sections in iV— space given in 
(12), the solutions NL[f a / P ](N, ji 2 ) can be used for further analysis directly. 

The next task is to compute the Mellin moments of the known coefficient functions for the 
different hard processes in case of massless fermions to 2-loop order which were usually computed 
in x— space. 

Let us start with the Drell-Yan process. We present here the relevant formulae for both 
unpolarized as well as polarized cross sections. The Drell-Yan process is given by 

H^Pi) + H 2 (P 2 ) - l+fa) + r(k 2 ) + X , (16) 

where Hi are the incoming hadrons with momenta Pj. fcj, with % — 1,2, are the momenta of 
final state leptons Z + Z~, respectively. Since the dominant contribution is through the 7-exchange 
s— channel processes 6 , we find the cross section can be expressed as 



d(A)a DY (x,Q 2 ) 4na 



2 



(A)W DY (x,Q 2 ), (17) 



dQ 2 3N C Q 2 S 
where (A) denotes the polarized case, Q 2 is the invariant mass of the di-lepton, 

s = (p 1 + p 2 ) 2 , x = %-> ( 18 ) 



3 The Z-exchange processes can be incorporated in a similar way. 
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N c is the number of colors, and a is the electromagnetic coupling constant. The scaling function, 
usually called the hadronic Drell-Yan structure functions (A) W DY (x, Q 2 ), is related to the 
coefficient functions as 

(A)W DY (x,Q 2 ) = ]T f dx,f dx 2 f dz{A)f a/Hl {x 1 ^ 2 ){A)f b/H2 {x 2 ^ 2 ) 



X 



(A) (z, QV) 5(x-zx lX2 ) , 



(19) 



where we have set the renormalization scale to be equal to the factorization scale. Alternatively, 
we can express the above equation as 



(A)W DY (x,Q 2 ) = £ f ld -^{A)^ ab {y^ 2 ){A)A^ [ X - )Q 2 ^ 

a,b= q ,q, 9 Jx y \y 

where the flux $ ab is defined by 

(A) $ ab (y, ^ 2 ) = I' - (A) f a/Hl (z, /i 2 ) (A) f b/H2 ( V -, p 
J y z \z 



(20) 



(21) 



The Mellin moment of or 1 \d (A) a DY (x, Q 2 )/dQ< 



reduces to 



M 



Id {A) a 



DY 



x dQ 2 



(A^ 2 ) = 



A-KQ 2 

3N C Q* 



M 



{A)W UY (iV,/i 2 ) , 



where 



M [(A) W DY ] (N, ^ 2 ) = M [(A) $ a6 ] (TV, /i 2 ) M [(A) A™] (N, /i 2 ) 



(22) 



(23) 



The flux M [(A) <E> ab ] (N,fx 2 ) can be computed using (15) with the input parton densities as 
boundary conditions. 

We now present the total cross section for both scalar (H) and pseudo-scalar (A) Higgs 
bosons at hadron colliders. The process is given by 

H 1 (P 1 ) + H 2 {P 2 ) - B + X , (24) 

where B = H,A. The total cross section for Higgs boson production is found to be 

7iG 2 B r 1 f 1 

^tot(^™ 2 ) = 8(N2 _ Yl J x dx ± J^ dx 2fa/ Hl (Xu ^)fb/H 2 (x 2 , ^) 



a,b=q,q,g 

xA ab R 



x 



2 ..2 



XiX 2 



-,m ,fi 



B — H,A, 



where m 2 is the mass of the Higgs boson and the scaling variable is defined by 



m 



x = ^-, S=(P 1 + P 2 ) 2 . 



(25) 



(26) 



A a b,B(x, m 2 , fJ 2 ) is the partonic coefficient function. Again we have set both renormalization and 
factorization scale to be equal. The overall constant G B is 



G B = -2 5 / 4 a s (/4) G 2 F r B F B (r B )C B (a s (/4), /4 m t 2 ) , 



(27) 



with 



a s (/4) 



a- 



Alt 



(28) 



where /ir is the renormalization scale. G F is the Fermi constant and the functions Fb(t) are 
given by 

F h (t) = l + (l-r)/(r), F A (r) = /(r)cot/3 

4m? 
r = — T 



/(r) = arcsin 2 — =, for r > 1, 



1 / i_vT^7 V 

— - I In , + m ] t < 1, 



(29) 



where /3 is the mixing angle in the two-Higgs doublet model. m t is the mass of the top quark. 
The coefficient Cb is given by 



C„ (a,(i/ R ), mf) = 1 + af(A) [bC A - 3C F ] + (afHp 



2 M 
2 F 



im C A C F + l ^C\ - ic F T f - \C A T S + (7C\ - UC A C F ) In ^ 



3 36 A 3 1 6 

+n ,7) ( - 4C F - ^CU + 8C F In 



m 2 



(a s (/i^),m 2 ) = 1, 
[26], where refers to the five-flavor number scheme. The color factors are 



(30) 



C A = N c , 



(31) 



Eq. (25) can be expressed in a compact form as 



1 B ( nG B l' ld y^ ( 2\V A 

) = ^3T) E X 7*<^ a; 

V / a,b=q,q,g Z> 



X 



I 11 
ub B I -• '/' •/' 



(32) 



The Mellin transform of (32) becomes 



M 



'tot 



TlG 



(N, m 2 ) = g(iy2 g ^ £ M [$ ab ] (iV + 1, ^ 2 ) M [A ab , B ] (iV, m 2 , //) 



ab=q,q,g 



(33) 



The flux M [$ a 5] can be extracted from the solution of the evolution equations (12) and M [A^s] 
can be computed analytically from the known functions A abjB (x, m 2 , /i 2 ), see appendix 3,4. 
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4 Finite Harmonic Sums 



We are dealing with the total scattering cross sections for the unpolarized and polarized Drell- 
Yan process and hadronic (pseudo) scalar Higgs boson production in the heavy-mass limit, which 
depend only on two variables, x = Q 2 /s and Q 2 , the invariant mass squared of the final state. The 
coefficient functions upto NNLO contain a large class of rr-space functions [7]. Their complexity 
reaches nearly 80, the number of all alternating and non-alternating finite harmonic sums [7, 
17, 18] of weight w < 4. The class contains simple functions of the form (1 ± re) -1 , ln m (x), 
ln m (l — x), and (ln m (l — x)/(l — x)) + ,m = . . . 3. More complicated examples are weighted 
Nielsen-integrals S njP (x)/(l f i). The Mellin transforms of these functions were evaluated in [7]. 
In the following we present some examples. 

The Mellin moment of the function 1/(1 + x) can be computed as follows: 



Jo 



dxx 



N-l 



1 + x) 



dxx 



f 

Jo 

(-If- 1 



N-l 



E(-i)v = E 



i=0 
■N-l 



(-1)' 

i + N 



i=i 



(-l) Ar - 1 [5_i(JV-l)+ln(2)](34) 



Likewise one obtains 
• l 

dxx 

o 



N-l 



(l-x) + 



dx 



x 



N-l 



- 1 



E 

i=0 



1 -x 

1 1 



I oo 

dx{x N ~ 1 - 1)E^ 



=0 



i + N i + 1 



-S 1 (N-1) 



where the finite harmonic sums are defined as 



A? 



N ( _ 1]k 



1=1 " 1=1 
Similarly, we compute the Mellin moment of x r ln(l — x) by expanding 



(35) 



(36) 



ln(l — x) 



St 



using partial fractions and changing the limits of the sums : 

r i °° 

/ dxx N - 1 x r \n(l-x) = -Y-r 
Jo y ' f^tiii 



E 



1i(i + N + r) N + rf^\i + N + r 



(37) 



N+r 

E 



S 1 (N + r) 



N + r — i 



N + r 



One may compute the Mellin moment of the distribution 



dxx 



N 



/ ln(l-x) ' 
I 1 -x 



dx 



x N ~ l - 1 



ln(l-x) N-l f 1 



1-x 



j: 



X 



N-2 



(38) 



ln 2 (l -x) 



N 



1 oo 

2 E-777 



iU^+j + JV-l) 



(39) 
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The summation is performed after partial fractioning and shifting the summation limits and one 
arrives at 



Jo 



dxx 



v-i ( H 1 ~ x i 
\ l-x 



= S^N-l) . 



(40) 



The latter sum is a multiple (nested) harmonic sum defined by 



N 



'm,tni,iri2 



,-,m k (N) — 



<7 

u mi, 



m2,...,m k V 



, Vmj > . 



(41) 



i=l 



In general, also alternating nested harmonic sums contribute, which are labeled with as well 
negative indices ki < 0, 



S, 



" [sign^)]"* ^ [sign(fc 2 )]"a 

2^ ifcii 2^ 



ni=l 



n2=l 



l*a| 



"jg 1 [sign(fc m )] f 

n m =l 



| &m | 



(42) 



with 7^ 0. Similarly one computes Mellin moments of more complicated function. For 
example, 



M 



dn 3 (l-x)' 
l-x 



(AO = lst(N - 1) + |s?(JV - 1)S 2 (N - 1) + ^Si(N - 1) 



+2S' 1 (JV - l)S- 3 (iV - 1) + -S 4 (N - 1) , 



(43) 



where the right hand side contains only single finite harmonic sums. The Mellin moments of 
Nielsen integrals can be mostly done by relating them to simpler moments after an integration 
by parts: 



Si dxx^S^-x) = 1| - ^ M[ ln 2 (l + x) (N) 



2N 



(44) 



The Mellin transforms of individual functions in rr-space may contain complicated sums [7], 
which cancel in combinations. As an example we mention 



dxx 



(-1) 



N-l 



Si, 2 (-x) + Li 2 (-x) ln(l +x) + - \n(x) ln 2 (l + x) 



N-l 



N 



s-ww + 1 ism - s-^N)]} + 1 - ± ^r 



N-l 



Cs ln(2)C 2 



(45) 



in which the sum S'_i i i(A r ) does not occur unlike the case for M[Li 2 (— x) ln(l + x)](N). We 
finally mention recursive integral-representations for the finite harmonic sums and weighted 
power sums, cf. [7], 



S± k (N) 



1 dx\ 

X\ 



(±x k ) N - 1 



dx\n k 1 (x) 



" (±x) k 



k=i 



k l 



(k-iy. Jo 

(l-l)\Jo 



(±x) 



N 



dzln 1 ' 1 ^) 



x =f 1 
±z) N - 1 



(46) 
(47) 
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The translation of the 2-loop Wilson coefficients [5,6] from x- to iV-space is performed using 
the tables given in [7] for the individual functions. One observes the cancellation of a series of 
sums. In particular we would like to mention, that multiple sums with the index { — 1} do not 
occur in the final results up to w=4. 

5 Algebraic Relations 

Finite harmonic sums obey algebraic equations, see e.g. [17]. The more simple relations for finite 
harmonic sums can be found in [7, 15,16,27]. Finite harmonic sums obey a shuffle algebra. To 
further simplify the expressions for the Wilson coefficients dealt with in the present paper the 
following relations are used for sums with two 

S ltl (N) = i(^(iV) 2 + S 2 (N)) 

5-1,1 (JV) = -S^iN) + SiWS-iiN) + S- 2 (N) 

S-i,- 2 (N) = -S^iN) + S^ 2 (N)S^(N) + S 3 (N) 

Si, 2 (N) = -S 2>1 (N) + S 2 {N)S 1 (N) + S 3 (N) 

S-i, 2 (N) = -5 2 ,_i(7V) + S 2 (N)S^(N) + S- 3 (N) 

51, - 2 ( N ) = S^ 1 (N) + S. 2 (N)S 1 (N) + S. 3 (N) 

S 2 , 2 (N) = ±(S 2 (N) 2 + S 4 (N)) 

S- 2 ,- 2 ( N ) = \(S^ 2 (N) 2 + S 4 (N)) 

5 2 , - 2 ( N ) = S- 2 , 2 (N) + S- 2 (N)S 2 (N) + S- 4 (N) 
S-i,- 3 ( N ) = -5 , -3,-i(^) + 5 , _ 3 (JV)S , -i(JV) + S , 4(JV) 

Si, 3 (N) = -S 3 , 1 (N) + S 3 (N)S 1 (N) + S 4 (N) 
S-i, 3 (N) = - ( 53,_i(7V)+ ( S 3 (7V)5_i(7V) + ( 5_ 4 (7V) 
S 1 ,_ 3 (N) = -S^ 1 (N) + S^(N)S 1 (N) + S^(N) (48) 

and three indices : 

Si,2,i(W) = -2S 2 , 1:1 (N) + S 3il (N) + S^S^N) + S 2:2 (N) 

Si,i,2{N) = S 2X1 (N) + [^i(7V)(5i !2 (7V) - S 2tl (N)) + S lt3 (N) - S 3>1 (N) 

Si,-2,i(N) = -2S- 2 , ltl {N) + S- 3A (N) + S 1 (N)S- 2 ,i(N) + 
Si,i,- 2 ( N ) = S- 2 ,i,i(N) + S- 2 (N)S 2 (N) - S- 2 , 2 (N) - S- 2 (N)S 1A (N) 
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+S 1 (N)S 1 ,- 2 (N) + S lt - 3 (N) - S 1 (N)S- 3 (N) 



(49) 



At intermediate stages of the computation we encounter various more complicated sums such as 
5*1,-1,2, S-i-i-2, S-i-2-i, S-2,-1,-1, 5*2,-1,1, 5i,2,-i, 52,1-1, 5_i,i,2, 5-1,1,2, 5_i,_3, and 5_i,3. All 
these sums do finally occur in symmetric combinations, which are polynomials of single harmonic 
sums. At the level of 3-fold sums only 5_ 2 ,i,i and 5 2 ,i,i remain. 



6 Coefficient Functions in AT-Space 



Let us now express the Wilson coefficients in terms of the remaining harmonic sums in terms of 
polynomials and as rational functions in A. Furthermore, we list the respective harmonic sums 
in terms of their Mellin-transforms, cf. [7] . 
The single harmonic sums are : 



5_ 4 (A) 

5_ 3 (A0 

5- 2 (A0 
5-i (AO 
5 4 (A0 

5 3 (A) 

5 2 (A) 

5i(A) 



(—1) N+1 - M 
6 



ln 3 (x) 



1 + x 



V ; 20 s ' 



(-lV v -M 

(_l)^+i M 
(-1) W M 



\n\x) 
1 + x 

\n(x) 
1+x 



(N + !)--& 



(N + !)--& 



1 + x 



i M 

6 



1 

— M 

2 



\n\x) 
l-x 

In 2 (a; 



(A + 1) - ln(2) 
2 



1 — x 



(A + 1) + jjCf 



(A+1) + Cs 



M 



ln(x) 
l-x 



-M 



1 — x 



(A+l) + C 2 



(A + l) . 



(50) 



These harmonic sums can be solely expressed in terms of Euler's ^-function and the /3- 
function [27] and their derivatives, which is related to the former combining two ^-functions 
with shifted argument. These functions represent at the same time the analytic continuation of 
these harmonic sums : 



S k (A) 



f— 1 ) k+1 

L_y_^(*-i) 



(N + l) + ((k) 



(k-l)\ 

S_i(A) = (-l)^(A+l)-ln(2) 



S- k (N) 



(fc-1)! 



pl*-»(N+l)- (l- — )Ck, k>2, 



(51) 
(52) 
(53) 
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(54) 



The following five double sums occur 

Li 3 (x) 
_ 1 + x 

1 



S_ 3A (N) = (-1)"M 



(N + 1) + ( 2 S_ 2 (N) - CsS-^N) - ^Cl + 2Li 4 (± 



+ ^3ln(2)-ic 2 ln 2 (2) + ^lu 4 (2) 



5_2,i(^V) = -(-l)^M 



S- 2 ,2( N ) = "(-1) JV M 



Li 2 (x) 
1 + x 



(N + l) + ( 2 S_ 1 (N)--( 3 + ( 2 \n(2) 



— |— (2Li 3 (x) - ln(x)(Li 2 (x) + C2)) 



(N + 1) + ( 2 S- 2 (N) 
1 



+2C 3 ^-i (iV) | ^ • ll.i, Q - ?Cs ln(2) + C 2 In 2 (2) ~ -h 



5 2 ,i(iV) = M 



Li 2 (x) 
1 -x 



(N + ^ + C^N) 



-Im 

2 



Li 2 (x) ln(x) 



1 — x 



Two triple sums contribute : 



S-2,i,i(N) = -(-1) JV M 



f 



1 + x 
1 



(N + l) + ( 2 S 2 (N) - ^(N) - ±S A (N) - ^C 2 2 (55) 



(N + 1) + CsS-^N) - Li 4 (i) + |cl + ^Csln(2) 



+ i C 2 ln 2 (2)--ln 4 (2) 



S 2A ,i(N) = M 



' Si, 2 (xY 

l-x 



(N + l) + C 3 Si(N) . 



(56) 



For the analytic continuation it is sufficient to determine M.[f(x)](N) for complex values of N 
since 

^ M[f(x)](N) = M [ln fc (x)/(x)] (N) (57) 

is easily obtained analytically. Therefore we will not count the associated derivatives as genuinely 
new functions. Besides Euler's ^-function the Mellin transforms of five further basic functions, 



LigOzQ Li 2 (^) S 1>2 (x) SiA x ) Li 3 (x) 
1 — x ' 1 + 1 — x ' 1+x' 1 + x' 



(58) 



are sufficient to express the different Wilson coefficients dealt with in the present paper. The 
Mellin transforms of these functions were calculated in [12] and are denoted by A k (N) with 
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k = 18, 3, 21, 8, 18, 6. In the appendices we furthermore refer to the functions A 5 (N) and A 2 2(N), 
which are given by 



A 22 (AO 



A 5 (N) 



_d_ 

dN 
_d_ 

dN 



A 3 (N) , 



A 18 (N) . 



(59) 



(60) 



The analytic continuation to iV e C is unique [28]. The hadronic cross sections are obtained 
calculating the inverse Mellin transform of (23,33) by a numerical contour integral 



Here the parameter c is the intersection of the contour and the real axis and is chosen right to the 
rightmost singularity of the function ~M[F(x)](N). The shape of the contour can be deformed 
at our convenience covering all singularities of M.[F(x)](N). 

7 Conclusion 

We have systematically analyzed the mathematical structure behind the NNLO coefficient func- 
tions for the unpolarized and polarized Drell-Yan process and hadronic scalar and pseudo-scalar 
Higgs boson production using Mellin moment techniques. Use of various algebraic and struc- 
tural identities, which relate the finite harmonic sums, reduces the complexity of the results from 
around 80 functions to only five basic functions, the -^-function and a few derivatives thereof. 
This is very useful both for the understanding of the nature of higher order corrections in the 
massless case and yields expressions which allow to perform fast numerical calculations at high 
precision for phenomenological applications and fits to data. The same structures are found in 
the case of polarized and unpolarized 2-loop fragmentation functions [21]. Together with the 
results of [20] for the Wilson coefficients for unpolarized and polarized deeply inelastic scatter- 
ing it is now shown that these structures are in common for all known massless 2-loop Wilson 
coefficients. 
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A.l Unpolarized Drell-Yan Coefficient Functions 

In this appendix we present the Mellin transforms of the coefficient functions for the unpolarized 
Drell-Yan Process. We used the same notation as in Refs. [5] and refer either to single harmonic 
sums Sk{N) or special functions, Ai(N). To avoid the occurrence of boundary terms in case of 
S^k(N) we use the functions P^ k ~^(N), k > 1, instead, with Af e C, outside the respective 
single poles of the functions below. 

A ab = Ai?+a s Ail ) +^Ai 2 ) (1) 



A (( 9 = 1 



(2) 



A^ — Cf 



(^-^- 8 ^-H in (K 



4 , R sm 



(ivTi7 + 16C2 - 16 + 8 (N + i) 



s ^ n + 1 KssKn-i) 



(3) 



A (1) = T f 
qg J 

4 



2 

+ T7 + 



6 



(N + l) N (N + 2)/ \fi 

4 



(91) jL 

: " ' 2 l + N 2 + (N + 1) ' N 



1 

+ TT 



Sm S 1 (N + 2) 
(N + l) 2 (N + 2) ' ~ (N + l) ' (N + 2) 2 N (N + 2) 



7 +8 S 1 (N+1) + 



(4) 



(2),S+V _ r 2 
*qq — 



18 + 32 S 2 (N - 1) + 32 S 2 (N - 1) - 32 ( 2 - 48 S ± (N - 1) 



x 



In 2 (Qpj + ^ (l28 - 192 S 2 (N - 1) + 64 ( 2 ^J S^N - 1) + 48 S 2 (N - 1) + 176 ( 3 



+48 S 2 (N - 1) - 93 + 24 C 2 - 64 Sf(N - 1) - 128 S 3 (N - 1)\ In 



256 C 3 



+256 S 3 (N - 1) )Si(N - 1) + - 64 ( 2 - 128 + 192 S 2 (N - 1) S 2 (N - 1) 



+32 Sf(N - 1) + ( - 64 C 2 - 128 ] S 2 (N -l)- 70( 2 - 60 Cs + ^ + 192 S 4 (N - 1) 



+ -C2 2 + 96S 2 2 (N-1) 



+CaCf 
-24Cs + 



f^-i)-n)m 2 g) + ((i6C 2 -f)^-i) 

193 88 



f S 2 (iV-l)-f 5 a (JV-l))ln(g) + (-56C 



176 



16 



+f S * iN - 1) + - 1) + ( - 166 + f - 1) + f Sf<* - 1) 



536\ 



12 



352 



| - I6C2 + — )S 2 (N- 1) --( 2 2 + 28( 3 + —S 3 (N-l) - — ■ — , 



1535 592 > 

+ ^c 2 



+C F rif 



2l'Q 2 \ , A6oa r ^ 34 , 16 o2 , 



2 _ ^(jv - 1) In 2 p- + — S 2 (N - 1) - — + — Sf(N - 1) 



-f C2 -^3(^l) + f-f^-l)-f^-l)-f^-l) 



^G + f Cs+(40A 18 (iV-l) + ^C 2 -16C3-^5 2 (iV-l) 



+16S 3 (iV-l) .Si(iV-l) 



52 446 



218 2278 
+ 



3N 2 27N 3(iV+l) 2 27(iV + l) 



-2AA 21 (N-1) + 



22 



22 



3(iV + l) 3N 



+ zrr= ln z ^ + 



Q 2 \ , I6C2 , 68C2 24C: 



+ 



N 2 3N N 



I6C2 20C 2 C 3 ^3(^+1) 4 ^(AQ 88£ 2 (iV+l) 

(jV+1) 2 3(iV+l) (iV + 1) (N + l) N 3(N+1) 

88 ,S'f( A' ) / , , C2 ... , -S' 2 i A* ) I 7G I .2 \ . S 2 (iV + 1) 



+ 



3iV V ^ ^ 3JV 2 9JVJ u ; \ 



(N + l) 



-48 ^ + 956 



(N + l) 9(N+1) 3(N+1) 



176 ^ Sl(iv + 1) + 16 M^ + 16 .-Miv + i) 



N 



(N + l) 



208 



+ - 28 C 2 + 28^(^-1) )Sf(N-l)+ - — + 4( 2 )S 2 (N-1) + 



8 



(iV+l)< 



32 



3(N+ 1) 



1 5 2 (JV + 1) + ( - A _ ^ 52 (7V) + 1 C2 2 + ((16 C 2 



-16 5 2 (7V-l))5 1 (-l + 7V) + 16C3-^-C 2 -8-^ T ^ ____ + i6_ 



496 88 140 o/Ar ^ 

+ ™ + ^ 2 (7V-i) + 



9(iV+l) 3N 2 3 



3(iV+l) 2 9N 



^-WA 18 (N-1) 



17 



88 S,(N) C2 
3N 



( W+1 ,-»t-»^-)>(s; 



, / 24 40 5i(iV + l) 8 



-:!2^4^-32 5 2 (AT-l) 



24 



, ,'g 2 \ 104 72 



iV 
48 



A^ 2 AT 



+ : ^^ T -216 1 S 4 (iV-l)-104 1 S 2 2 (iV-l)- + 



48 



(iV + 1) 2 (N + l) 
-16A 21 (N-1)- ^ 



100 / C2 
+ 232 



iV 3 (iY+1) 3 
112 



(iV+1) (iV + 1) 2 



+ 



192 



(iV+1) 3 
C2 



104 



A^ 4 (N + l) 4 
S2(N + 1) 12 \ ( , C2 24 C 2 24 C 3 



C2 24 C 3 , 24 ^3^+1) , 21 ^ 3 (iV) s A 18 (iV) 



-120 — - • Hi , 

(AT + 1) 2 (N + l) (N + l) (N + l) 



N 



N 



A 1; (N + 1) t 64 192 _ W & _ M?U (JV) + ( 288 & 

(iV + 1) \ iV 2 iV 3 AT JV JV y \ 



-288 S 2 (iV - 1)) S 2 (iV - 1) + 64 + 6 4 ■+ 



156 



A" 



(N + l) 



(N + iy 



+ 



64 



(N+l) 



)sKN+l) + f^+[-^-^) S!(N) + ((-256 C 2 



+256S 2 ( J V-1))S 1 ( J V-1) + 



144 If \„„, „ ,144 



+ iH) SlW + 3 2^M+^i 
N N (N + l) 



, 48S 2 (AT-l) + 32 



ff 2 (N+l) 56 
(A" + 1) N 2 



72 (T 72 

+32A 18 (N-l)-—+A8C 2 + imS,(N-l)-9Q^- w —^ 



96 



-160Ca-96 



+ 



16 



ff 2 (iV+l) 
(N + l) ' (N + l) ' N 



120 96 ^ V 



8OC2 S 2 (AT-1) + - 



16 



+ 



60 



(JV + 1) (AT + 1) ; 



N 
112 + 



(TV) 



+ 1^256 ( 3 — 256 5* 3 (A^ — 1) — 32 A l8 (N — l)j Si(N — 1) 



18 



A _ A -r — C put 

qq,AA t J 

8 



2 



In 2 ^ + 



+ 



3N 3{N+\)) \fi 
16 



16 ^(7V + 1) 32 \ /Q 2N 
9AT ' 9(JV+1) ' 3 AT 3(iV+l) 2 3 (iV+1) 3^/ ^// 2 , 



16^(7V) 



/ 16 32 \ /Ar . 40 8 , 
8C2 16 S?(N) 



4 188 20 40 
3i^ + 27iV" SiV 3 ^ 3 



+ ; 44 ^ - - A 18 (7V - 1) + - 5 3 (AT - 1) - - 20 XQ 
3N 3N 3(N+1) 2 3 8V ; 3 3V ; 3(N + l) 3 



16 Sf(N + 1) 



8C2 



412 



3(N + 1) 3(N + 1) 27(N+1) 3(N + 1) 



8 g 2 (AT + l) 8 S 2 (AT) 



3N 



A {2) - = A (2) - 

qq,AA qq,AD 



CAC F -^f) 



32 94 
-24C2 + 18C3- + "T7 + 



30 



78 



•18 S 3 (iV - 1) - 2 S 4 (iV - 1) + 6 Si(N -!)+_ + 



N 2 AT (N + l) 2 (N + l) 
8 



N 3 (N + l) 3 



12A 18 (N-1) 



12 



+32 A 21 (N- I) + — + 



12 



/ 16 16 

, + 12A 22 (N- 1) + 

N 4 N + l 4 22V J \ N 3 N 2 



16C 2 + 16S 2 (N-1) U 2 (N-1) 



C2 , -, „ C3 - C2 , -, r. C3 



+ 4 



S 2 (N + 1) 
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+ 4 - 12 *£± 1} - 12 ^ - 8 A )^ N t 1] S^ + ((-l Q S 2 (N 



N 



N 



(N + l) 

+ 16(2)3^-1) - 16A 18 (N- 1) + 12(2 + 



„ S 2 (N+1) g g 2 (N) | 
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28 



AT 



8 +8^1 



(N + l) (N + l) 3 N 



32 



(iV + 1) 



at 1 (nTW + Ni + ^-u^n-i)- AT 



(ivTT) ) ln (f ) " X ^ + ( 24 Ca + 12 * 2(iV " 1} " 24 Ss{N - 1} 



19 



-12C 2 )si(iV-l) + 



16 



56 S 2 (N + 1) g C2 



(AT + 1)3 (N + l) (N + l) 



(N + l) 



16 



\S 1 (N + 1) 



(10) 



A (2) _ r 
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16 



40 

+ + 



+ 



16 



3N 3 N 2 3N 3(AT+1) 3 3(N + l) 2 3(N + 2) 3 
40 



(AT + 2) 2 3(N + 2) 



A (2) _ = A W _ 



32 (3^(N + 1) 64 /3W(iV + 2) 32 /3«(3 + AT) 



3N 



3(JV + 1) 
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36 54 
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26 24 
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32 



+ 
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+ 48 
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A. 2 Polarized Drell-Yan Coefficient Functions 



In this appendix we present the Mellin transforms of the coefficient functions for of the polarized 
Drell-Yan process. We used the same notation as in Refs. [5]. 
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A. 3. Scalar Higgs Coefficient Functions 

In this appendix we present the Mellin transforms for the coefficient functions of the hadronic 
Higgs-pro duct ion cross section in the heavy-mass limit. We used the same notation as in 
Refs. [6]. 
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A. 4 Pseudo-scalar Higgs Coefficient Functions 

In this appendix we present the Mellin transforms of the coefficient functions of the hadronic 
pseudo-scalar Higgs-boson cross section in the heavy-mass limit. We used the same notation 
as in Refs. [6]. 
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